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Abstract 

We discuss an interacting tachyonic dark energy model in the context of the holographic principle. 
The potential of the holographic tachyon field in interaction with dark matter is constructed. The 
model results are compared with CMB shift parameter, baryonic acoustic oscilations, lookback 
time and the Constitution supernovae sample. The coupling constant of the model is compatible 
with zero, but dark energy is not given by a cosmological constant. 
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In the last years, there have been several papers where an interaction in the dark sector 
of the universe is considered A motivation to considering the interaction is that 

dark energy and dark matter will evolve coupled to each other, alleviating the coincidence 
problem A further motivation is that, assuming dark energy to be a field, it would 
be more natural that it couples with the remaining fields of the theory, in particular with 
dark matter, as it is quite a general fact that different fields generally couple. In other 
words, it is reasonable to assume that there is no symmetry preventing such a coupling 
between dark energy and dark matter fields. Using a combination of several observational 
datasets, as supernovae data, CMB shift parameter, BAO, etc., it has been found that 
the coupling constant is small but non vanishing within at least la confidence level 
U Q], in two recent works, the effect of an interaction between dark energy and 
dark matter on the dynamics of galaxy clusters was investigated through the Layser-Irvine 
equation, the relativistic equivalent of virial theorem Using galaxy cluster data, it has 
been shown that a non vanishing interaction is preferred to describe the data within several 
standard deviations [10]. However, in most of these papers, the interaction term in the 
equation of motion is derived from phenomenological arguments. It is interesting to obtain 
the interaction term from a field theory. Some works have already taken a step in such a 
direction jg], 11]. On the other hand, there have been several papers where the dark energy 



is associated with the tachyon scalar field. The tachyon field has been studied in recent 
years in the context of string theory, as a low energy effective theory of D-branes and open 



strings 



12| . The pressure of the tachyon fluid is negative, and it has been used in cosmology 



as a candidate to dark energy 8], 13] - 16]. The first question about tachyons concerns 
the choice of the potential. Common choices for the tachyonic potential are the power law 
and the exponential potentials, both capable of reproducing the recent period of accelerated 
expansion, the last of these being motivated by some string theoretical models. However, 
these choices are in fact arbitrary. In principle, any other form for the potential which leads 
to recent accelerated expansion would be acceptable. 

On the other hand, it is possible that a complete understanding of the nature of dark 
energy will only be possible within a quantum gravity theory context. Although results 
for quantum gravity are still missing, or at least premature, it is possible to introduce, 
phenomenologically, some of its principles in a model of dark energy. Recently, a combination 
of the tachyon model with the holographic dark energy model has been made available [if]] - 



previously, combinations of quintessence and quintom models with holographic dark energy 
had been proposed [nj], jl8]. Specifically, by imposing that the energy density of the tachyon 
fluid must match the holographic dark energy density, namely pa — 3c 2 MpjL -2 , where c is 
a numerical constant and L is the infrared cutoff, it was demonstrated that the equation 
of motion of tachyons for the non-interacting case reproduces the equation of motion for 
holographic dark energy. In fact, to impose that the energy density of tachyons must match 
the holographic dark energy density corresponds to specify the potential of tachyons. This 
can be seen as a physical criterion to choose the potential. Here, we generalize this idea for 
the interacting case. 

We consider the lagrangian density 

Cm (x) = V=g~ [-V(<p)y/i - ad^d^ + '-[^YV - $ V m7 ^] - (M - /fy>)M j , 

where a is a constant with dimension MeV~ 4 , /3 is the dimensionless coupling constant, 
V(f>) the potential and g the determinant of the metric. From a variational principle, we 
obtain 

17*%* - M*# = , (1) 
z(V M *) 7 M + M** = , (2) 



where M* = M — (3cp, and 



gMV^g)gV 1 dlnV{<p) 
V M cr<£i + a— — — h 



(3^ 



a/1 - a&hpd^p . 



(3) 



- ad^ipd^ip a dip aV(<p) 

The eq. (0Q) and (J2]) are, respectively, the covariant Dirac equation and its adjoint, 
in the case of a non vanishing interaction between the Dirac field and the tachyon 
field ip. For homogeneous fields and adopting the Friedmann-Robertson- Walker metric, 
g^=diag(l, — a 2 (t) , — a 2 (t) , — a 2 (£)), where a 2 (t) is the scale factor, the eq. (QQ) and (J2]) 
lead to 

d(a 3 ^ty 



dt 







and (131) reduces to 



oap 



a 



oap 1 



dip aV(<p) 
where H = - is the Hubble parameter. From (J3J) we have = ^o^o (~)' 



(4) 
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From the energy-momentum tensor, we get 

V(cp) 



^ a/1 — oap 2 
P v = -V(<p)y/l-a<p 2 

p* = o . 



(6) 
(7) 



Notice that from (GO) and (JEJ) we have u; ¥ 
to time and using (jSJ) and (HI), we get 



3i 



a<p 2 — 1 . Deriving (jBJ) and (j7|) with respect 



p v + 3H Pv (ojp + 1) = /3<p^o^c 



aoy 
a / 



and 



p^ + 3Hp^ = — ^(p^o^o 



oo\ ! 
a / 



(8) 



(9) 



where the dot represents derivative with respect to time. The Friedmann equation for a flat 
universe reads 



H 2 



1/2 

where Mpi = (8iiG) is the reduced Planck mass. 

In order to determine the dynamics of the interacting tachyon, it is necessary to specify 
the potential V(ip). In [8], a power law potential had been chosen. Here, instead of choosing 
an explicit form for V(ip), we will specify it implicitly, by imposing that the energy density 
of the tachyon fluid, given by ([6]), must match the holographic dark energy density, p\ = 
3c 2 M 2 Pl L- 2 , where c is a numerical constant and L is the infrared cutoff. The evolution 
of the interacting tachyon fluid with redshift will be given by the equation of evolution for 
the holographic dark energy density, with a certain expression for the equation of state 
parameter u v . In fact, we will see that imposing the energy density of tachyons to match 
the holographic dark energy density leads to an expression for the potential of tachyons. 



3M* Pl 



a 



+ 



V(<p) 



a<p z 



(10) 



In 



19| it has been argued that, in order that holographic dark energy drives the recent 



period of accelerated expansion, the IR cutoff L must be the event horizon Rh. Substituting 



Rh in the expression of the holographic dark energy, we get Rh 



dt' 



therefore, 



a it') a it) ' 



Differentiating both sides with respect to time, using the Friedmann equation (|T0|) together 
with conservation equations (JHJ) and (jUJ), we obtain 



dQ v 
dz 



2a, 



l + z 



p 'iVl^piOp + 1 



(11) 



Equation ffTTj) is just the equation of evolution for the holographic dark energy [19 |. 



We define r = — . Deriving r with respect to time, using (jSJ) 



and eliminating (p by 



using 



±./ 1+UJ * we obtain 



r = 3HruJ v — sign [ip] 



/3(l + r) 2 v /TT 



(12) 



3Mj >l ^H 2 VI + z o, 

The sign of ip is arbitrary, as it can be modified by redefinitions of the field, ip — > —ip, and 
of the coupling constant, — > — /3. We can rewrite ^o^o in terms of observable quantities. 
In fact, by imposing that the dark matter density today matches the observed value, we 
obtain M^o^o = — f °p — M 00 , where we defined = i/^V 9 - Furthermore, noticing that 



00 



after some algebra 



we can eliminate r and r in favor of O^ and O v in (fl2|) . Using (TITI) we obtain, 



w^, (z) 



1 

3 



2^/0777) 7 ( z ) 
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where 



with 



E(z) = 



H(z) 



7W + 
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z) — 0o and 5 = — — is an effective coupling constant. Notice that, 



where A0 (z) = 

if 5 = 0, ( fT3|) reproduces the equation of state parameter obtained in 19]. 
The evolution of the tachyon scalar field is given by 

d(j) _ a/1 + ^4, (z) 



(15) 



dz H E (z) (1 + z) 

From pip and f|T5|) we can calculate the evolution with redshift of all observables in the 

model. If we wish to calculate the time dependence, we need to integrate the Friedmann 

equation ( TTOl) . which can be written in the form 

dt 1 

~d~z ~ ~ H E(z)(l + z) ' 



Here it is worth saying that in the holographic dark energy model, in the non interactin 



case - f|T3|) with 5 = - u$ can be less than —1. However, as already mentioned in 16], 
if we wish that the holographic dark energy is the tachyon, then because ( II 5p . must be 
more than —1. Nevertheless, in the interacting case considered here, due to the fact that 
depends explicitly on 0, can not be less than —1. On the other hand, the square root in 



( IT3|) must be real. We can verify that ou<p is real and > — 1 if (i) — — < 1 or (ii) ^ — > 1 

and > 2 1 n ^ - ^ _ ^ _ However, case (ii) is irrelevant, as it corresponds to large 

values of -j|L For example, if f^o = 0.7 and c = 0.8, we have > 2.45. Below, we will 
see that the observational data constrain ~ 10 -1 . In order that be real for all future 
times, as Q,^ — > 1, it is necessary that c > 1. 

It is interesting to notice that the condition "^ 0O < 1 is precisely the condition for which 
the entropy of the universe increases jl9]. As fi^, — > 1 in the future, it is necessary that 
c > 1. Therefore, the condition for be real is precisely the same one for the entropy to 
increase. So the model respects the second law of thermodynamics. 

From ([6]), we can compute the potential V (z) as 



V(z) 



e (z) a 



(Z) yJ-U^z) , (16) 



PcO 

where p c o = SM^Hq, E (z) is given by ( 1T4|) . (z) is given by ( |T3l) and (z) is the solution 
of ffTTj) . From ([TBI) and f|T5|) . we can compute V(<ft). In figured! V(<f)) is shown for some 
values of 5 and c. Notice that, as we chose positive, then evolves to the mininum of the 
potential. However, if we had chosen <p negative, then because the right hand side of f|T5|) 
would has the opposite sign, V(<fr) would be now an increasing function of and again 
would evolve to the mininum of the potential. 

The equation for evolution of (|T5l) can be written in an integral form as 



Since the model depends on A0 - through E (z) - and neither on nor on 0o, then it is 
independent of 0o- In other words, 0o is not a parameter of the model and can be chosen 
arbitrarily. Therefore, the parameters of the model are 6, c, h and Q^o- Below, we discuss 
the comparison with observational data and the results obtained. 



In 



20] . the lookback time method has been discussed. Given an object % at redshift Zi, 



its age t(zi) is defined as the difference between the age of the universe at Zi and the age of 




FIG. 1: Potential of tachyons V ((/)), in units of p C Q = SM^Hq. <f> is in units of H^ 1 . The solid 
lines are for c = 0.85, the dashed ones are for c = 1.1 and the dotted are for c = 1.35. For each 
value of c the curves from right to left are for = —0.1 (red), = (black) and = +0.1 
(green), respectively. 

the universe at the formation redshift of the object, zp, that is, 

dz' r°° dz' 



t(zi) = H, 1 



[l + z')E(z>) 



z F 



;i + z')e(z') 



dz' 



(l + z')E(z 
where ti is the lookback time, given by 



- = t L {Zp) - t L {Zi) , 



, S rr 1 f" dz' 

Using (IT7j) . the observational lookback time t°£ s {zi) is 

= - = [*f - *(*)] - K bs - t L (z F )} 

= C - t(Zi) - df , 

where iff 5 is the estimated age of the universe today and df is the delay factor, 

i.obs 



We now minimize Xiw 



df = t° bs -t L (z F ) 



i=l 



(17) 



(18) 



where ^(z^p) is the theoretical value of the lookback time in Zi, p denotes the theoretical 
parameters, t°^ s (zi) is the corresponding observational value given by (j!8j) . Uj is the uncer- 
tainty in the estimated age t(zj) of the object at Zi, which appears in (lisp and a t ou is the 
uncertainty in getting t^ s . The delay factor df appears because of our ignorance about 
the redshift formation zp of the object and has to be adjusted. Note, however, that the 
theoretical lookback time does not depend on this parameter, and we can marginalize over 
it. 

In 21( and 22j the ages of 35 and 32 red galaxies are respectively given. For the age 



231 ] . Although this estimate for t 



obs 



of the universe one can adopt tg bs = 13.73 ± 0.12Gyr 
has been obtained assuming a ACDM universe, it does not introduce systematical errors in 
the calculation: any systematical error eventually introduced here would be compensated 
by the adjust of df, in ffl8|) . On the other hand, such an estimate is in perfect agreement 
with other estimates, which are independent of the cosmological model, as for example 



t. 



obs 



12.6t.2^Gyr, obtained from globular cluster ages 24j and t 



obs 




12.5 ± 3.0Gyr, 



obtained from radioisotopes studies [25]. 

For the cosmic radiation shift parameter in the flat universe we have 

■ dz' 



R= 



M 



E(z' 



where z\ s is the last scatte ring surface redshift parameter. The value R has been estimated 



231 ] results as R Q b s = 1.715 ± 0.021, for the flat universe, with 



from the 5-years WMAP 
Zi s = 1090.5 ± 1.0 and is very weakly model dependent 



261 ] . Thus we add to \ 2 the term 



XCMB 



(R-R, 



obs) 



a 



R 



Baryonic Acoustic Oscilations (BAO) are described in terms of the parameter 

2/3 



-1/3 



1 


f ZBAO dz' 1 


_ZBAO J 


o E(z')\ 



where zbao = 0.35. It has been estimated that A b s 
X 2 the term 



0.493 ± 0.017 \M. We thus add to 



Xbao 



(A -A 



obs) 



<J, 



The BAO distance ratio r BAO = D v (z = 0.35) /D v (z = 0.20) = 1.812±0.060, estimated 



from the joint analysis of the 2dFGRS and SDSS data 



281 ] . has also been included. It was 
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demonstrated in [28j that this quantity is weakly model dependent. The quantity Dy (zbao) 
is given by 



D v {zbao) = c 
So we have the contribution 



z BAO 



zbao dz t 



X 



H(z BAO )\Jo H(z') 
_ (r-BAo - rflo) 2 



1 1/3 



rBAO 



a 



rBAO 



Finally, we add the 397 supernovae data from Constitution compilation 29J. Defining 
the distance modulus 

" z dz' 



H{z) = 5log 10 
we have the contribution 



cl+2 



E{z')\ 



+ 25 — 5logioH 



397 



2 _ ^ Z i) ~~ Vobs{Zj)] 2 

xsn — 2_^i 



i=i 

Using the expression x 2 = Xm + Xcmb + Xbao + xL AO + Xsn> the likelihood function is 



given by 



C(S, c, h, Qfo) oc exp[- 



X 2 (S, c, hfClfa) 



In table 1 we present the values of the individual best fit parameters, with respective la, 
2a and 3a confidence intervals. 

Table 1: Values of the model parameters from lookback time, BAO, CMB and SNe la. 



.5 

H Q 


-0.066 ± 0.092 ± 0.185 ± 0.277 


C 


D oce:+0.056+0.141+0.260 
U.OOO_ .022-0.035-0.045 




0.720 ±0.016 ±0.033 ±0.049 


h 


0.666 ±0.013 ±0.026 ±0.038 



The figures [2] and [3] show the marginalized probability distributions for 5 and c. The 
coupling constant 5 is compatible with zero within la level. Therefore, in this work it 
was not found evidence for interaction. However, the best fit value is of same order of 
magnitude obtained in a previous work, where the interacting tachyon model with power 
law potential has been considered [81]. f^o corresponds to a density matter parameter 
today Q m0 = 0.280 ± 0.016, in perfect agreement with cosmological model independent 



estimatives, as for example flMobs = 0.28 ± 0.04 [30j. The value of h = 0.666 ± 0.013 is 
also in excellent agreement with observational values, independent of cosmological model 



{h obs = 0.69 ± 0.12 [2l| and h obs = 0.72 ± 0.08 |3l|). 




5/H 



FIG. 2: Probability distribuction of 5. 




n — 1 — i — 1 — r 

0,90 0.95 1,00 1,05 1.10 1.15 1.20 



FIG. 3: Probability distribuction of c. 



Figure H] shows the joint confidence regions for two parameters. As we can see, there is 
little degeneracy between the parameters of the model. In the confidence regions for 5 versus 
c and for c versus fl^o, we see that there is a lower limit on c ~ 0.8. This also can be seen 
in the marginalized probability distribution of c, which dies for c < 0.85. This lower limit 



is explained by the condition 



< 1, necessary for to be real and > —1, discussed 
10 
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— i — i — | — i — | — i — | — i — | — i — | — i — | — i — | — i H 1 1 1 1 1 1 — i — i — i 1 1 1 1 1 — 

-0,4 -0,3 -0,2 -0,1 0,0 0,1 0,2 0,3 -0,4 -0,3 -0,2 -0,1 0,0 0,1 0,2 0,3 



8/H 5/H 




FIG. 4: Confidence regions of la, 2a and 3a for two parameters. 

above. This limit can be seen more clearly in c versus f^o confidence regions. Moreover, 
we have c ~ y^TJio for the best fit values of these parameters. This implies that u^o ~ — 1 
and the model approaches ACDM today. This is consistent with the fact that, as ACDM 
fits all observational data, then any alternative model must not deviates much from ACDM 
for z^0. However, for z > 0, the model is qualitatively different from ACDM, as 
approximates —1/3, see figure [5) 

We have obtained c < 1 at la confidence level. As already said above, this implies 
that the equation of state parameter will not be real for all future times. However, 
this is not a very serious problem, because c is compatible with values above unit at 
2cr confidence level. Moreover, one could say that c < 1 is only an effect due to lack 
of more precise observational data. Anyway, the very simple model presented here is 
expected to be only an effective description of a more sophisticated subjacent theory of dark 
energy. In principle, nothing guarantees that it will be a good description for all future times. 
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-0,3- 
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-0,5- 
-0,6- 
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ln(a) 

FIG. 5: Equation of state parameter of dark energy w^(a), for c = 0.85 and = —0.1 (red dashed 
line), 77^ = (black solid line) and = +0.1 (green dotted line). 

In summary, a combination of holographic dark energy model and interacting tachyon 
field was implemented. It was showed that it is possible to fix the potential of interacting 
tachyon by imposing that the energy density of the tachyons must match the energy density 
of the holographic dark energy. A comparison of the model with recent observational data 
was made and the coupling is consistent with zero. However, in a previous work [8| with 
interacting tachyons with power law potential, a non-vanishing coupling constant had been 
obtained with 90% of confidence. So the possibility of a small, but calculable, interaction 
in the dark sector remains open, and future investigations, with more realistic models and 
more observational data are necessary to solve this question. 
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